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We study properties of some known trial wavefunctions in bilayer quantum Hall systems at the 
total filling factor vt = 1. In particular, we find that the properties of a meron wavefunction and 
a natural "quasi-hole" wave function are dramatically difi'erent due to the broken symmetry and 
the associated Goldstone mode in the bulk. Although the ( smallest ) meron has localized charge 
1/2 and logarithmically divergent energy, the charge of the quasi-hole excitation extends over the 
whole system and its energy diverges linearly with the area of the system. This indicates that 
the natural quasi-hole wavefunction is not a good trial wavefunction for excitations. It also shows 
that the energy of the naive candidate for a pair of meron wavefunction written down previously 
increases quadratically instead of logarithmically as their separation increases. Our results indicate 
that qualitatively good trial wave functions for the ground state and the excitations of the interlayer 
coherent bilayer quantum Hall system at finite d are still not available and searching for them 
remains an important open problem. 

I. INTRODUCTION 

The variational approach has been successful for investigating condensed matter systems, especially in supercon- 
ductivity and Quantum Hall Effects. In 1957, Bardenn, Cooper and Schrieffer proposed BCS trial wave-function for 
a 5'-wave superconducting ground state and also constructed all possible quasi-particle excitations above the ground 
state. In 1983, Laughlin proposed a trial-wave function for the ground states of Fractional Quantum Hall Effects 
(FQHE) [1] at filling factors v — -^j^- He also constructed trial wavefunctions of quasi-hole and quasi-electron 
excitations and found that they are localized and carry fractional charges zfc 2s+i ■ J^-i^ proposed Composite Fermion 
(CF) trial wavefunctions for the ground states at Jain's principal series at filling factors v — 2sp±i PI- 
ory produces a spectrum of the low energy CF excitations which matches that obtained from exact diagonalization 
study in finite size systems. These excitations are localized and carry fractional charges 2sp±i ■ Halperin generalized 
Laughlin's wavefunctions to {m,m',n) wavefunctions ( See Eqn.22 ) for FQHE in a two-component system [3]. These 
components could be the spins of electrons when the Zeeman coupling is very small or layer indices in multi-layered 
system. Quasi-hole and quasi-particle wavefunctions for generic (m,m' ,n) state in two-component systems can also 
be constructed by analogy to those in the single layer systems ( see Eqn. 23). Scarola and Jain constructed a general 
class of CF wavefunctions in bilayer Quantum Hall systems which include Halperin's (m, m', n) states as a subset and 
investigated the energies of these wavefunction as a function of the distance d between the two layers [4] . 

One of the most interesting cases is the spin-polarized Bilayer Quantum Hall systems (BLQH) at total filling factor 
vt = 1. For sufficiently small interlayer separation d, Halperin's (111) wavefunction ( sec Eqn. 6 ) is believed to be a 
good trial wave-function for the ground state. In contrast to general {m,m' ,n) state, the (111) state ( or {m,m,m) 
state in general ) is a symmetry broken state and there is an associated gapless Goldstone mode in the bulk [5-7]. 
1^ • This system has been a subject of enormous experimental and theoretical investigation over the last decade [8-10]. 
^ I When the interlayer separation d is sufficiently large, the bilayer system decouples into two separate compressible 
O • J/ = 1/2 layers. However, when d is smaller than a critical distance dc, in the absence of interlayer tunneling, the 
] system undergoes a quantum phase transition into a novel spontaneously interlayer coherent incompressible phase 
^ [10]. Recently, at low temperature, with extremely small interlayer tunneling amplitude, Spielman efc al discovered a 
k>( \ very pronounced narrow zero bias peak in this interlayer coherent incompressible state [11]. M. Kellogg et al also 
d ' observed quantized Hall drag resistance at h/e^ [12]. 

I By treating the layer index as a pseudo-spin index, Girvin, Macdonald and collaborators mapped the bilayer system 
into a Easy Plane Quantum Ferromagnet (EPQFM) [10] and explored many rich and interesting physical phenomena 
in this system. The low energy physics is given by an effective 2 + 1 dimensional XY model. They interpreted 
the gapless mode as the spin wave mode of quantum XY model, and the vortex excitations of the XY model as 
topological defects called " merons " which carry fractional charges ±1/2 and also have ± vorticities. The merons 
have logarithmically divergent self energies and are bound in pairs at low temperature. The lowest energy excitations 
carry charge ±e which are a meron pair with opposite vorticity. By drawing analogy with superconductivity, they 
wrote down the EPQFM wavefunction Eqn. 7 which can also be interpreted as the pairing of an electron in one layer 
and a hole in another ( exciton condensation ) ( see also Ref. [5] ). They also wrote down the wavefunctions for the 
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( smallest ) mcron Eqn.l4. Recently, Ye applied both Mutual Composite Fermion ( MCF) and Composite Boson 
(CB) Chern-Simons theories to study imbalanced i^t = ^ bilayer [13,14] and trilayer quantum Hall systems [15]. He 
found that CB approach is more superior than the presently known MCF approach in the study of interlayer coherent 
multi-layer QH systems. 

Because there is a gap in the bulk for generic {m,m! ,n) states, both the ground state wavefunctions ( See Eqn.22 
) and the quasi-hole wavefunctions ( see Eqn. 23 ) are good trial wavefunctions. However, in the interlayer coherent 
Quantum Hall (m, m, m) states, there is a neutral gapless mode even in the bulk. Earlier exact diagonalization studies 
[16] showed that (111) wavefunction ( See Eqn.22 ) remains a good one in the d ^ limit. However, as pointed out 
in [10], the (111) wavefunction may not be qualitatively good at finite d, because (111) is a broken symmetry state 
of isotropic ferromagnet at d = instead of a easy-plane ferromagnet at finite d. Unfortunately, a qualitatively good 
ground state wavefunction at finite d is still unknown [17]. In this paper, instead of finding such a good ground state 
wavefunction at finite d, we focus on the quasi-hole wave functions ( see Eqn. 23 ) built on (111) state and measure 
its energy above the (111) state. So far their properties have not been investigated in detail. Their relation to meron 
wavefunctions in second quantization form is not clear. For example, are they good trial wavefunctions ? Are they 
new excitations different than merons ? The relation between (111) wavefunction and the EPQFM wavefunction in 
second quantization form is also not explicitly established. 

The gapless mode is a hallmark of the interlayer coherent Quantum Hall state. Its existence may dramatically alter 
the properties of the wavefunctions of the ground state, quasi-hole and quasi-particle. In this paper, we will study 
the relation between (111) wavefunction ( Eqn. 6 ) and the EPQFM wavefunction ( Eqn. 7 ). We find that they are 
equivalent. We will also investigate the properties of meron wavefunction ( Eqn. 14 ) and the quasi-hole wave function 
( see Eqn. 12 ). These two wavefunctions differ only by "normalization factors" ( the "normalization factor" is defined 
just above Eqn. 2 ). Normalizations factors have been shown not to be important in single layer Quantum Hall ( 
SLQH ) systems ( see section HA). However, we find that they make a dramatic difference in BLQH. Although the 
smallest meron has a localized charge 1/2 and logarithmically divergent energy, the charge of the quasi-hole excitation 
extends over the whole system and its energy also diverges linearly as the area of the system size. This indicates that 
the quasi-hole wavefunction is not a good trial wavefunction for any low energy excitations. We also find the energy 
of the possible wavefunction of a pair of merons ( Eqn.32 ) written down in [10] increases quadratically instead of 
logarithmically as the separation of the pair increases. 

The rest of the paper is organized in the following. In Sec. II, we will discuss the relations between the (111) 
wavefunction and the EPQFM wavefunction, the quasi-hole wavefunction and the meron wavefunction. In Sec. HI, 
we will study the charge distributions and energies of the meron and the quasihole, first by plasma analogy, then 
by Monte-Carlo simulations in both disk and spherical geometries. We will also discuss the properties of the trial 
wavefunction of a pair of merons. Finally, we state the conclusions in Sec. IV. 



II. WAVEFUNCTIONS IN BILAYER QUANTUM HALL SYSTEMS 

We will first review some subtleties in the trial wavefunctions in single layer QH system, then point out these 
subtleties make important difference in the charge distributions and energies of different trial wavefunctions in BLQH. 



A. Subtleties in the trial wavefunctions in single layer QH system: 

In symmetric gauge, the single-particle eigenf unctions of kinetic energy and angular momentum in the LLL are 
At filling factor u = l/m where m = 2s -|- 1 is an odd integer, the Laughlin ground state wavefunction is: 

N 

j<k ki,---,kN 

where X^^i ki = ^^-^^^^^—^m which is the total angular momentum of the state. 

If we insert one infinitesimally thin flux quantum at position z = 0, then the single particle state is changed from 
(f>m{z) to (/)m+i{z) = ^2(^ _)_^) '^'"(-^)- wli^re the prefactor will be called normalization factor in the following. 

The Quasi-hole wavefunction can be written as [18] ( In the following, we will suppress the Gaussian factor ): 
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V2^(A;i + !)••• (few + 1) 

Removing all the normalization factors in front of the coefficients Cfei,...,fc„, the above wavefunction is simplified to 

fci,-",A;jv i 

which is the well known Laughlin quasi- hole wavefunction. It carries localized fractional charge l/rn. 

Note that for Integral Quantum Hall Effect (IQHE) m = 1 when tpi is a single Slater determinant, it can be shown 
that all these normalization factors factor out of the determinant, therefore, removing these factors has no effect at 
all, namely, V'm.g/s ^^'^ i^m,qh are identical when m = 1. However, for FQHE m > 1, V'm.q/i ^^'^ '4'm,qh are different. 
Obviously, both and tpm.qh carry the same fractional charge 2m+i ■ SLQH, it was shown by numerical studies 

in finite size that tpm,qh and have similar energies [19]. However, tpm,qh reduces the density near the origin 

without ruining the good correlations in the ground state wavefunction tpm, its energy is expected to be slightly lower 
than that of V;^,,^. 

'4^m,qh can be best derived from Jain's CF theory of SLQH system where the ground state wavefunction Eqn.l is 
rewritten as [2]. 

N 

■^m{Zl, - ■ ■ ,Zn) = Y[{Zj - Zkf^Xl{zi,---,ZN) (4) 

where the 2s fold zeros {zj — Zk)"^^ incorporates the strong correlations of electrons in the ground state and xi is just 
the wavefunction of one filled Landau level. 
The wavefunction of a quasihole is: 

N 

■4'm,qh = Yl i^J ~ Zkf^Xl,hole, Xl,hole = (H 
j<k i 

where Xi,hoie is the wavefunction of a hole at = 1. 

If the QP and QH are far apart, the wavefunction for the QH at the origin is exactly 4'm,qh instead of V'm,?/!' because 
inserting one flux quantum affects only xi without touching the correlation factor [zj — Zfe)^^. 

Unfortunately, except for IQHE m = 1, there are no simple second quantization forms for ipmj^m,qh,4'm qh- 
conclude that in SLQH system, the energies of trial wavefunctions are insensitive to whether or not we keep track 
of the normalization factors due to the existence of the gap in the bulk. In the following, we will apply the insights 
gained from SLQH systems to BLQH systems and arrive at very different conclusions. 

B. Ill state in first and second quantization: 

In first quantization, the ground state trial wavefunction of BLQH in d — > limit is Halperin's (111) state [3]: 

Ni N2 Ni N2 

where z and w are the coordinates in layer 1 and layer 2 respectively. 

In the EPQFM picture [10,5], the ground state EPQFM wavefunction was written in second quantization form: 

|G>= n ^{ci,+ci,)\Q> (7) 

m=0 ^'^ 

where M = N is the angular momentum quantum number corresponding to the edge of the system. 

The wavefunction Eqn.7 can also be interpreted as the pairing of an electron in one layer and a hole in another 
which leads to exciton condensation. In the following, we will show that the EPQFM wavefunction in the second 
quantization form of Eqn.7 is equivalent to (111) wavefunction in Eqn.6. 
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We can expand Eqn.7 into: 

1 ^ 1 

JVi=0 ^' fei,...,fejv 

where N = Ni + N2 and P is the permutation of N variables. 

Projecting onto the state with A^i electrons in layer 1 and N2 electrons in layer 2 ( this state will be denoted by 
{Ni,N2) in the following ) and then transforming into the first quantization form leads to: 

<f>k.,+Awi) ■■■^k. {WN,){1, T) • • • (iVi, T)(iVi + 1, i) • • • {N, i)] (9) 

where A stands for anti-symmetrization. 

Moving the summation over the orbital states into A, the sum is essentially the (111) state: 

^ (^^) . . . ^^^^ {ZN,)ct>k^^ + Awi) ■ ■ ■ <t>k^{wN,) 

Ni N2 Ni N2 

= n n^^' ~ '^j)!!'^^' ~ ^jOH'^"'' ~ "^j) = ipni{z,w) (10) 

i=lj=l i<j i<j 

Finally 

\G; NuN2 >= {l=)^^^-^A[Mz,w){^, T) • • • (iVi, T)(iVi + 1, i) • • • {N, i)] (11) 

This concludes that the orbital part of the projection of the EPQFM wavefunction Eqn.7 onto (A^i, A^2) sector is 
exactly the same as (111) state. In the following, we will study the excitations of the (111) state. 



C. " Quasi-hole " versus Meron wavefunction: 

In first quantization, the most obvious trial wavefunctions for the " Quasi-hole" at the origin in layer 1 and layer 2 
are [10]: 

^gh,iiz,w) = (Y[zi)^ui, ^qhAz^w) = (YIuh)-^ in (12) 

i=l i=l 

We can get the above trial wavefunctions by naively generalizing the composite fermion quasi-hole wavefunction 
Eqn.5 in SLQH to the corresponding quasi-hole wavefunction in BLQH: 

Zi - Wj)Y[{wi - Wj) X (one hole in layer 1) (13) 

i=lj = l i<j 

In second quantization, the meron wavefunction in EPQFM picture [10], created by inserting one flux quantum in 
layer 2 at the origin, is 

M-l 

\Meron,2>= -T^iC^ + CU,0\^ > (14) 

m=0 ^ 

where 0, J, state is not occupied. 

An important question to ask is whether the second quantization form of the meron wavefunction is equivalent to 
the "quasi-hole" wavefunction in the first quantization form in Eqn.l2. 

In the following, we will show that they are different. 

We can expand Eqn.l4 into: 
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1 ^ 1 

V ^ ;vi=0 

Projecting onto the state {Ni, N2) and then transforming the above equation into the first quantization form lead 
to: 

|meron,2;7Vi,7V2> =(i=)^^^^ J2 (-ir('=--''=-)-^^[</.fe, (^1) • • • </)fc,, (^ivj 

fci ,---,kN 

<^fc„,+i+iK) • • • <^fc„+i(«;ivj(l, T) • • • (iVi, r)(iVi + 1, i) • • • {N, i)] (16) 
Using (pm+iiw) = I ^ (t>m{w) and moving the summation over the orbital states into A, we find the orbital 

y 2(m+l) 

part becomes 

X =(t>kM, + liWl) 7^5=TT<^fc«("^^2) (17) 

-y/2(fcjvi+i + 1) \/2{kN + 1) 

Naively, one may expect that all the normalization factors can be factored out just like in m = 1 IQHE in SLQH. 
However, detailed calculations show that this is not the case. In fact, these normalization factors play an important 
role. 

If we neglect the normalization factors just like what we did in SLQH in Eqn.3, the above sum leads simply to the 
quasi-hole wavefunction in layer 2 listed in Eqn.l2. 

N2 

'^qh,2{z,w) = {Ylwi)'^lll{z,w) (18) 
i=l 

If keeping all the normalization factors, then we note that 4>m{z) peaks strongly around its maximum at \z\^ = 2m. 
If \wi\ > 1, we can replace . in Eqn.l7 by t^. 

Therefore, far away from the origin, we have: 

^meronA^, W ^ Oo) = (V] -^)*iii(2, w) (19) 

Let us look at the following complete wavefunction which includes both spin part and orbital part 

-( ' 
where Z = {z,w). 

Expanding and then projecting onto (A^i, A^2) sector and extracting the orbital part, we can see the orbital part is 
exactly the same as Eqn.19. 

Setting Zi = |Zi|e*^', the spin part of the wavefunction can be written as [10]: 

XW = ^(>) (21) 

which is clearly a vortex. 

It was shown that the meron wavefunction Eqn.l4 or 17 has logarithmically divergent energy with the system size 
Em ~ InM where M = N is the angular momentum quantum number corresponding to the edge of the system [10]. 

This shows that the orbital part of the projection of the second quantized form of meron wavefunction Eqn.l4 
onto {Ni,N2) sector is not the same as the quasi-hole wavefunction Eqn.l2. The properties of the "quasi-hole " 
wavefunction Eqn.l2 have not been seriously studied so far. In the following, we will study the charge distributions 
and energy of this wavefunction and find that they are quite different from those of meron wavef unctions. 



*=nif z./k>^(^) (20) 
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III. INVESTIGATION OF THE ENERGY OF QUASI-HOLE WAVEFUNCTION 



In this section, from plasma analogy, we will study the charge distributions and energies of quasi- hole wavefunctions 
qualitatively and gain some physical insights. Then we will stiidy them quantitativ(^ly by Monte-Carlo simulations in 
both disk and spherical geometries. We will also compare with the meron wavefunctions. 

A. Plasma analogy of the quasi-hole wavefunction 

We will first review Halperin's {m,m',n) state [3] and then study the very interesting case (111) state which is a 
interlayer coherent state in this series. The (m,m',n) wave function is given by: 

Ni N2 Ni N2 

i=l j=l i<j i<j 

The quasi-hole wavefunctions were written down in [10] : 

'bgh,l{z,w) = (Y[zi)'^m,m',n, *g/i,2(-2, w) = (]^ Wi)*m,m',n (23) 

i=l i=l 

These trial wave functions clearly reduce the density near the origin without ruining the good correlations in the 
ground state (m, m',n). The fractional charges of these quasi-hole states were studied in [10] from two component 
plasma analogy. From the perfect screening conditions, the fractional charges of ^qh,i{z,w) in layer 1 ( ) and in 
layer 2 ( 63 ) are found to be: 

„* _ * _ "'^ *_ — n m' + n 

7 2' ^2 — '1 2' — '/ 2' — '/ 2 \-^^) 



mm — n 

where , e* are the total and relative charges. 

It is easy to see el and 63 always carry opposite charges as expected because of the strong interlayer correlations 
built in the (m,m',n) ground state wavefunction. It is also easy to see = i^i where i/i is the filling factor in layer 
1. The fractional charges for 'i'qh,2{z, w) can be obtained from the above results by m ^ m', v\ — > 1^2 where V2 is the 
filling factor in layer 2. 

For the (111) ground state, it is easy to see both ej and 63 diverge, while = vi remains finite to be equal to the 
filling factor in layer 1. In the balanced case z/i = 1^2 = 1/2, el — 1/2, namely, the total charge is 1/2 [20]. However, 
e* diverges. It is well known that there is a capacitive energy term Ec ^ J cP"r{p2 — piY for BLQH with finite distance 
d. This term strongly suppress the relative density fluctuations between the two layers. Therefore, the divergence of 
e* indicates that the energy of quasi-hole may diverge much faster than that of meron wavefunction which diverges 
only logarithmically [21]. 

In the following, we will study the charge distributions and energy quantitatively by Monte-Carlo simulations in 
both disk and spherical geometry. 

B. Monte-Carlo simulations 

The charge density operator is defined as p{z) = Yld=i ^{^ ~ ^i)- SLQH system, its average value at z is 

p{z) = N J d''z2 ■ ■ ■ d^ZN\^{z, Z2,---, zn)\VZ, (25) 

where Z = J d'^z\ip{z)\'^ and the normalization is / d^zp{z) = N. 

For BLQH system, the charge densities at the two layers pi{z),p2{w) are: 

pi{z)=Ni j d^z2---d^ZN, j d''w\^l){z,w)\^/Z 

P2{w) = ^2 j d^z j d'^W2 ■ ■ ■ d'^WN^ \il){z, w)f/Z (26) 
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with Z = J zd?"w\ij){z,w)\'^ and the normaUzations J (Pzpi{z) = Ni,J (Pwp2{w) = N2. 
In the foUowing, we focus on the balanced case Ni = N2 = N/2. 
For the ground state (111) wavefunction, we find: 



^ M-l 

PIO(^) = P20{Z) = Po{z) = 2 Z] (27) 





While, for the meron wavefunction, we find: 

M-l , M 



Plm{z) = \'t>m{z)\'^ = Po{z), p2m{z) = ^ \^m{z)f (28) 

^ ^1 

The deviations from the ground state in the two layers are 

Spim{z) = plm{z) - Po{z) = 0, Sp2m{z) = p2m{z) - Po{z) = -^(1^510(2:)^ " \(I>m{z)\'^) (29) 

The charges in top layer remains untouched, while 1/2 charge is moved from the origin to the boundary ( Fig.l ), 
namely, there is a hole with charge 1/2 at the origin and extra charge -1/2 at the boundary. 

I ' layer 1 

^ ^ A layer 2 

charge 1/2 

Fig 1: The charge densities of the smallest meron wavefunction 

The fact that the charges in top layer remains unmodified when one flux quantum is inserted in the bottom layer 
is because the meron wavefunction ignores the strong interlayer correlation. We expect the neutral gapless mode will 
strongly modify the charge distributions of a meron. 

We can compare the above meron wavefunction to the wavefunction of the ordinary charge 1 excitation: 

Ni N2 

'^chl{z, w) = 

i=l i=l 

Obviously, it has finite energy. Its charge distribution is depicted in Fig. 2. 

A charge 1/2 a 
^ layer 1 



f 



A I A 

charge 1/2 



layer 2 



Fig 2: The charge densities of the charge 1 excitation 
The density profile for the quasi- hole in layer 2 is: 

SP.M = ^ E™.,...,.„(m...i + l)-(m. + l) 

'''^^'^^^ = Y E..,...,™.(-...r + l)-(m. + l) 



(31) 



The above expressions can be evaluated analytically only for small number of particles. Fortunately, Monte-Carlo 
simulations can produce various quantities for large number of particles. In the following, we will first discuss the 
results in disk geometry, then in spherical geometry. 
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1. Results in disk geometry 




Fig 3: The charge densities of the quasi-hole wavefunction in layer 1 and 2 

As shown in Fig. 3, inserting one flux quantum at bottom layer pushes electrons towards the boundary to create 
charge deficits near the origin at the bottom layer. At the same time, electrons are pulled from the boundary to 
create charge surplus near the origin on the top layer. However, the excitations are not localized but extends over the 
whole system. Both density plots seem to collapse roughly on a single curve (except around the boundary) when the 
distance is scaled by the system size R = \/2N/ v, namely, Spiqh{r) ~ fi{r / R), 6p2qh{r) — f2{r/R). We can see at 
r/R ^ 0.64 both density deviations change sign. This value is close to the radius r/R = ~ 0.7 where the area 

inside is equal to the area outside in a disk geometry. We may tentatively define this radius as the size over which 
the quasi-particle charge extends. 

We can evaluate the total "quasi-hole" charge pqht — ^Piqh + ^P2qh and charge density difference profile Spgh = 
Spiqh — Sp2qh between the two layers. They are shown in Fig. 4a and Fig. 4b respectively. Both extend over the 
whole system. Again Spqh = fi{r/R) — f2{r/R) = f{r/R) changes sign at r/R ^ 0.64. The fact that all pght curves 
approach the horizontal axis as N increases is expected, because the total area enclosed by these curves should be 
the total charge —1/2 which spreads over the whole system, therefore the total charge density approaches zero as N 
increases. 




Fig 4: The total and relative charge densities of the quasi-hole wavefunction 

The integrated total and relative charges inside the radius r, i.e. qi{r) = Jr,^rd^r'[pi{r') — po{r')],i = 1,2 are 
shown in Fig.Sa and 5b. 
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(b) 



A -0.3 







:;; - 




\ 1 - 


N= 50 




- 100 




200 




400 





0.2 0.4 0.6 0.8 
r/R 



10 I — 

— 

-10 - 

i-20 - 

r-30 - 

-40 - 



= 50 
100 
200 
400 



0.2 0.4 0.6 
r/R 



Fig 5: The integrated total and relative charges inside radius r. 

As shown in Fig. 5, the total integrated charges at r/R ^ 0.64 are indeed -1/2, which is consistent with the result 
achieved in the plasma analogy. However they do not show a plateau because there is no plateau in the relative density 
( Fig. 3) either. The integrated relative charges simply diverge linearly as the system size grows. This is expected, 
because Q{r) = qi{r) - q2{r) = 27r ,f{r'/Rydr' = f^'"' f{x)xdx = R^Fqir/R). Setting ro/R - 0.64, wc get 
Q{ro) ~ ii^ ~ A^. The plasma analogy implies the total charge difference Q diverges, the Monte-Carlo simulations 
explicitly show that Q diverges linearly as N. We can also estimate the energy of this quasi-hole by first estimating 

q2 

the capacitive energy Ec = where C is the capacitance which is inverse proportional to the distance d between 
the two layers and proportional to the area the charge spreads over. From the above analysis Q ^ N,C ^ N, then 
Ec ^ N. We expect the energy of the quasi-hole should be at the same order as Ec, therefore Egt ~ N, namely, the 
quasi-hole energy diverges linearly as the system size grows. This linear divergence is much faster than the logarithmic 
divergence of the meron energy. 

We also show the quantitative capacitive energy Ec = J d'^f{p2 — pi)^ in Fig. 6 




Fig 6: The capacitive energy of BLQH 

In Fig, 6 is shown that Edro/R) ~ 0.015A'', namely, increase linearly as A^. This is expected, because Ec{r) = 

2-K f{r'/R)r'dr' = 2-kR^ j^'"- f{x)xdx = R^Fc{r/R), Setting ro/i? ~ 0.64, we get Ec{ro) ^ R^ N which 
is consistent with the above estimate of Ec- The plasma analogy implies the quasi-hole energy Eqh diverges, the 
Monte-Carlo simulations explicitly show that Eqh diverges linearly as A^. 

Note that Eqn.30 is essentially the wavefunction of a charge 1 excitation in single layer spin polarized = 1 
Quantum Hall state. It can also be considered as a smallest skyrmion. When we split the skyrmion, namely, replacing 
the prcfactor (IlilJi Y[f=i ''^'i) iii Eqn.30 by Hil^il^i ~ ^o) nfc?i(^i ~ '"'o) and then gradually increasing |zo — wo\, 
its energy Echi ~ jzo — w^oP increases quadratically. Since, in the |2;o — m'oI — > cxd limit, Eqn.30 will reduce to the two 
charge 1/2 quasi-hole wave functions in Eqn.l2. 

An approximate wavefunction for a pair of meron which includes both spin and orbital part was proposed in Eqn.208 
of [10]: 
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where Z = (z, w). 

As discussed in [10], the spin is up at wq and down at zq. By symmetry, it seems that there may be 1/2 charge located 
at Wo and zq, so it may be a good candidate for the wavefunction of a pair of merons. By a unitary transformation ( 
Eqn. (110) in Ref. [10] ), Eqn.32 can be mapped to a skymion wavefunction with a finite radius. This shows that at 
zero distance, the candidate wavefunction for a pair of merons Eqn.32 has the same energy as that of a skymion and 
one meron can be viewed as half of a skymion. However, at finite distance, the Hamiltonian has only U{1) symmetry, 
the unitary transformation does not respect this f/(l) symmetry anymore. Eqn.32 has much lower energy than that 
of a skymion, since it keeps the spins in the XY plane. 

Now let's look at the ener^ of Eqn.32. It is easy to see that the projection of this wavefunction to a given sector 
(7Vi,7V2) is ni=i('2« ~ ^o)ni=i(^j — wo)^'iii, so its energy Epair ^ |zo — w^ol^ instead of logarithmically as naively 
expected, because the charges are extended between zq and wq. A good wavefunction for a pair of meron whose 
energy increases logarithmically as the distance between the two merons is still unknown. 

We can also calculate directly the energy of the (111) state Eq and the quasihole state Eq^. The intra- and inter-layer 
Coulomb interactions are given by: 

Vn = V22 = 1/r, V12 = V21 = 1/ + d? (33) 

where d is the distance between the layers. 

Unfortunately, the boundary contributions from the positive background charges in the disk geometry are not easily 
taken care of. In order to avoid this difficulty, we perform the simulations in spherical geometry. 

Through the use of the single-particle eigenfunctions (pq.m oc uQ+^yi-^e^Q'!' in the lowest Landau level, we can 
construct Halperin's (111) wavefunction in spherical geometry : 

Ni N2 Ni N2 

^niiz,w) = HHiiu^ij - ViU,)e'^^'+^^')] Y[[{u^VJ - v,Uj)e'^^^+^^^] UHu^d, - 5,S.)e*('^'+^^)] (34) 

j=lj=l i<j i<j 

where z = {u,v) with u = cos |e~"^/^, t; = sin|e"'^/^ are the spinor coordinates in sphere 1 and w = {u,v) with 

u = cos le"*"^/^, V = sin le*"^/^ are the spinor coordinates in sphere 2 [22]. 
The quasi-hole wavefunction at north pole in sphere 1 is 

^chAniz, w) = ([[ Vie"^-/^)^iii (35) 

i=l 

which corresponds to increase the monopole strength Q in the center of the sphere 1 by 1/2, namely, Q ^ Q + 1/2. 

Wc evaluate the energy difference Eh — Eq in spherical geometry where the background positive charge distributions 
cancel. The results at several finite separations are shown in Fig. 7. 
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Fig 7: The energies of quasi-hole wavefunction at diflferent separations d > of the two layers in spherical geometry 



As shown in Fig. 7a, Eh — Eq = C{d)N with C{d) ^ d ^ as d ^ 0. Wc also simulate the best fit for the power 
law behavior E^ — Eq = C{d)N^'^ at a given d as the system size N increases. As shown in Fig. 7b, the power Pd{N) 
are slightly larger than one in finite size, but all approach to one in the N ^ oo thermodynamic limit. 

We also calculate — Eq at d = 0. The result is shown in Fig. 8. 




Fig 8: The energy of the quasi-hole wavefunction at d = in spherical geometry 

In sharp contrast to all the results at > shown in Fig. 7, the energy at = remains finite in the thermodynamic 
limit, positive constant. This behavior is expected, because d = limit is a singular limit where the capacitive energy 
Ec vanish. 

We also perform the simulations of charge density distributions and capacitive energy in spherical geometry. The 
results are similar to those shown in Figs. 3-6 in disk geometry. 

From the results in both disk and spherical geometries, we conclude that the quasi-hole energy diverges linearly as 
the system size grows which is much higher than that of meron. 



IV. CONCLUSION 



In this paper, we investigated the properties of some known wavefunctions and their relations. The existence of the 
gapless mode in the interlayer coherent Quantum Hall state dramatically change the properties of the wavefunctions 
of the ground state, quasi-hole and quasi-particle. We studied the relation between (111) wavefunction ( Eqn.6 ) and 
the EPQFM wavefunction ( Eqn.7 ) and found that they are equivalent. We also investigated the properties of meron 
wavefunction ( Eqn.l4 ) and the quasi-hole wave function ( see Eqn.l2 ). These two wavefunctions differ only in 
normalization factors. Although normalizations factors are presumably not important in SLQH systems ( see section 
IIA) because of the gap in the bulk, we find that the charge distributions and energies of trial wavefunctions in BLQH 
are very sensitive to the normalization factors. This is because there is a neutral gapless mode even in the bulk in 
BLQH. We find that although the smallest meron has a well localized charge 1/2 and logarithmically divergent energy, 
the charge of the quasi-hole excitation extends over the whole system and its energy diverges linearly with N. This 
is because the gapless mode dramatically suppresses the relative density fluctuations between the two layers at any 
finite separations. Those results indicate that the quasi-hole wavefunction is not a good trial wavefunction for the 
excitations of (111) state. 

As pointed out previously in Ref. [9], (111) wavefunction may not be a qualitatively good wave function due to 
the existence of the gapless mode of the easy-plane ferromagnet at finite d. Unfortunately, the qualitatively good 
ground state wavefunction is still unknown. It is still unknown if the energy of the quasi-hole wave function Eqn.l2 
will become logarithmically if we replace the (111) wavefunction by the unknown qualitatively good ground state 
wavefunction. 

An important open question is to determine what is the lowest energy excitation: a pair of meron or conventional 
charge 1 excitation in Eqn.30 ? As shown in the text, the energy of the naive candidate for the wavefunction of a 
pair of meron Eqn.32 increases quadratically instead of logarithmically. Again, it is not known if the energy of the 
pair will become logarithmically if we replace the (111) wavefunction by the unknown qualitatively good ground state 
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wavcfunction. Without knowing a qualitatively good wavefunction for a pair of merons, it is hard to answer this open 
question precisely. 

In general, trial wave function approach is a very robust approach to study SLQH and multi-components systems as 

long as there is a gap in the bulk. The gap protects the properties of the system such as charge density distributions and 
energies from being sensitive to some subtle details of wavefunctions. However, the results of this paper indicate that 
qualitatively good trial wave functions for the interlayer coherent bilayer quantum Hall system at finite d, especially 
for the excitations, arc still unknown and searching for them remains an important open problem. Fortunately, 
effective theory approaches such as EPQFM approach [8-10] and Composite Boson theory approach [7,9,14,15,13] 
are very effective to bring out most of the interesting phenomena in this system. In fact, all these effective theories 
start from the insights gained from Halpcrin's (111) wavcfunction. Although many interesting results can be achieved 
qualitatively from effective theories without knowing the precise wavefunctions for the ground state and excitations, 
some quantitative results such as the lowest energy gap, quantum phase transition at dc, finite temperature phase 
transition, etc can be better achieved if we have a better knowledge of the microscopic wavefunctions. 

This work is partially supported by NSF under grants Nos. DMR-0240458 ( GSJ ). We thank J. K. Jain and A. H. 
Macdonald for many insightful discussions and communications throughout this work. 
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